New contribution to WIMP-nucleus scattering 



Gary Prezeau, Andriy Kurylov, Marc Kamionkowski, Petr Vogel 
Division of Physics, Mathematics, and Astronomy, 
California Institute of Technology, Pasadena, CA 91125 

A weakly-interacting massive particle (WIMP) is perhaps the most promising candidate for the 
dark matter in the Galactic halo. The WIMP detection rate in laboratory searches is fixed by the 
cross section for elastic WIMP-nucleus scattering. Here we calculate the contribution to this cross 
section from two-nucleon currents from pion exchange in the nucleus, and show that it may in some 
cases be comparable to the one-nucleon current that has been considered in prior work, and per- 
haps help resolve the discrepancies between the various direct dark-matter search experiments. We 
provide simple expressions that allow these new contributions to be included in current calculations. 

PACS numbers: 95.35.+d 



A neutral, stable, weakly-interacting massive particle 
(WIMP) with mass near the electroweak scale is one of 
the most natural dark-matter candidates P, because 
it has a cosmological density comparable to that con- 
tributed by halo dark matter. The most widely stud- 
ied WIMP candidate is the neutralino — a linear combi- 
nation of the supersymmetric partners of the photon, 
Z and Higgs bosons — in the minimal supersymmetric 
extension of the standard model (MSSM). There are, 
however, other possibilities, including the sneutrino (the 
neutrino's superpartner), heavy neutrinos, Kaluza-Klein 
modes in models with universal extra dimensions, etc. 

If such WIMPs exist, one way to observe them is di- 
rectly via detection of the O(30 keV) recoil energy im- 
parted to a nucleus in a low-background detector when 
a WIMP collides with the nucleus 0, Q. The pre- 
dicted event rate for these searches depends primarily 
on the cross section for elastic WIMP-nucleus scattering. 
Typically, a WIMP-nucleon interaction is derived from 
the WIMP-quark interaction that appears in the MSSM 
(or other WIMP-theory) Lagrangian. It has long been 
known that for supersymmetric models — and recently 
shown for more general pointlike WIMPs Q — the WIMP 
can couple to the nucleon with either a spin-independent 
(SI) or spin-dependent (SD) interaction. Nuclear physics 
then allows the WIMP-nucleus cross section to be derived 
from the SI and SD WIMP-nucleon interactions. 

In this work, we use an effective field-theory descrip- 
tion of nuclear interactions to consider the contribution 
to WIMP-nucleus scattering from two-nucleon currents, 
heuristically, the coupling of WIMPs to the virtual pions 
that hold the nucleus together. Although this contribu- 
tion is, as one might expect, small in many cases, for 
SI interactions it is usually non-negligible and sometimes 
comparable to or even bigger than the one-nucleon cur- 
rent usually considered. 

We begin with the underlying WIMP-quark interaction 
Lagrangian, construct effective hadronic operators, and 
organize them according to the order in qj Ah with Ah ~ 
1 GeV a hadronic scale like the nucleon mass or the chiral 
symmetry breaking scale and q ~ , the pion mass with 
q being the momentum transfer or the Fermi momentum 
of the nucleons inside the nucleus. We then focus on the 



neutralino to perform an explicit evaluation of the two- 
nucleon scattering amplitude and we estimate the range 
for its magnitude. 

The two-nucleon contributions will have a long range 
piece through the exchange of pions and a contact piece 
as shown in the diagrams of Fig.Q]where the WIMP is de- 
noted by X . The long range contributions are enhanced 
compared to the two-nucleon contact interaction as can 
be seen from power counting in the small momentum 
q. To show this, we first organize the A-hadron vertices 
in powers of q using chiral perturbation theory [(| and 
the fact that the momentum transfered by the WIMP is 
- 100 MeV or less for WIMP masses ~ 100 GeV or less. 
The leading order (LO) quark operators should therefore 
induce effective WIMP-hadron vertices that do not in- 
volve derivatives of the pion fields or pion mass insertions, 
(which are always quadratic and therefore do not con- 
tribute at LO or next-to-leading order (NLO)) the NLO 
operators would involve a single derivative of the pion 
field, the next-to-next-to-leading order (NNLO) would 
involve two derivatives or pion mass insertions and so- 
on Q. Thus, counting pion propagators as 0(1/ q 2 ) and 
the pion- nucleon vertex as 0(q), we find that the ma- 
trix element of Fig. ma) ~ 0(K wn q- 2 ), Fig.[Hb,c) ~ 
0(K NNn q- r ), Fig.Hd) ~ O{K m q ), where the Ki 
denote the order of the A"-hadron vertices. In gen- 
eral, the LO A"-hadron vertex in each diagram will have 
Ki <~ O(q ), though in certain cases symmetry consider- 
ations of the underlying particle-physics model require 
that the LO vertex vanish 0. Thus, the long range 
two-nucleon contributions of Figs. ^ a), and I2b,c) are 
enhanced by 1/q 2 and l/q, respectively, relative to the 
short-range operator of Fig-Ed)- In what follows, we 
will at most consider NLO contributions to the scattering 
amplitude and will henceforth neglect the contact term. 
Therefore, we consider all terms in and Knn^ to 
0(q) and 0(1) respectively. 

We now construct the A'-hadron effective Lagrangian 
from an underlying A"-quark effective Lagrangian where 
A" is a general neutral WIMP field with arbitrary spin, 
isospin, and parity (note that we are not yet special- 
izing to the neutralino). The most general set of CP- 
conserving interactions between WIMP and quark cur- 
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FIG. 1: Two-nucleon diagrams that contribute to WIMP-nucleus scattering where the WIMP is generally denoted by X. Graph 
(a) is of 0(l/q 2 ), graphs (b) and (c) are of 0(l/q) while the contact term of graph (d) is of 0(1). The exchange diagrams are 
not included. The filled circles represent the non-standard model vertices. 



rents is 



-Xq 



G F ^2 [a\Sqq + a q 2 Vq^q 
i 



+a q T^qa^q],(l) 

where the S, V, V M , A^ 1 and T^ v are respectively lin- 
ear combinations of scalar, pseudoscalar, vector, axial- 
vector, and tensor operators constructed from A and ar- 
bitrary constants. The constants a\ are parametrized in 
terms of the Fermi constant Gf — v / 2ff 2 /(8A/^) (where 
Mw is the VT^-boson mass) and depend on the particu- 
lar particle-physics model (such as the MSSM) used to 
generate the A-quark Lagrangian. The inclusion of n- 
body quark operators with n > 1 would induce further 
suppression in a heavy scale such as Mw and are ex- 
cluded. In the following, we will neglect the a\ term 
which is generally suppressed by an inverse power of the 
WIMP mass. Also, only SU(2) isospin symmetry will be 
considered since kaon-exchange currents are suppressed 
in elastic WIMP-nucleus scattering. We will therefore 
truncate the sum in Eq. to up and down quarks when 
considering vertices that contribute to exchange currents. 

The Lagrangian of Eq. (JIJ generates various A-hadron 
interactions vertices. For the particular problem of A- 
nucleus scattering at NLO, we will focus on the following 
vertices: ir 2 X 2 , ttN 2 X 2 (which contribute to the two- 
nucleon current at NLO) and N 2 X 2 (the one-nucleon 
contribution to the scattering amplitude). The corre- 
sponding A-hadron Lagrangian will have the general 
form 



c 



ttN" 



c 



(2) 



the last of these is that which has been considered 
in prior work, while the first two are new. The three 
parts of Cx in Eq. @ are characterized by parameters 
bi,Cj,dk, respectively, which in turn depend on af of 
Eq. Q). The precise shape of this general Lagrangian 
will therefore depend on the A-quark Lagrangian and 
we now proceed to construct each term in Eq. @ in turn. 

ir 2 X 2 vertex. Up to NLO in powers of the pion 
momentum, the ir 2 X 2 Lagrangian will look like 

E™ = b s n ■ ttS + b v i(w x d M 7f) 3 V„. (3) 

In Eq. ©, only the terms with two WIMPs will be kept 



in the expansion of S,V, V M , and A^ . To derive an ex- 
pression for b s , we use the matrix element, 



b s (TT a \Tf ■ Tf\TT a ) = G F (n a \ a,lqq\w a ), 



from which we obtain using soft-pion techniques 



6 S = 



GFml K + af), 



2(m u + m d ) 



(4) 



(5) 



where the m; are current quark masses. The sum of the 
quark masses is given in Ref. Q as 5 MeV < m u + md < 
11 MeV and we will use the average value of 8 ± 3 MeV. 

For the NLO term of Eq. , we use the conservation 
of the vector current (CVC) to write 

b v = G F (a u 3 -a d 3 ). (6) 

Note that in the MSSM there is no contribution from a*. 

ttN 2 X 2 vertex. The NLO two-nucleon currents will 
only receive contributions from the a\ 5 terms in Eq. 
(the other possible terms contribute at NNLO) yielding 



c v N^i{r x n) 3 NV„ + Ca N^j 5 i(f x tt) 3 AU m ,(7) 



where N is the nucleon isospinor. To extract the hadronic 
scales that appear as one matches the quark-A La- 
grangian to the hadron-A Lagrangian we can use dimen- 
sional analysis and the scaling rule 0, El 



NN \ k ( \ 1 ( 7T \ m G F 



An/^y ^Ah/ Viw a h 



x (A H /.) 2 , (8) 



where /„• = 92.4 MeV is the pion decay constant and 
(k, I, m) refer to the form of the hadronic part. The Cj's 
of Eq. J7J which have one nucleon current and one pion 
[hence, (fc,Z,m) = (1,0, 1)] can be rewritten 



^~*~F ( u d \ r 

c v,a = —j. (0^5 — a 4,5/"v,a i 



(9) 



where S v is O(l) and from CVC we have S a — — ,9a/2 
with <7a = 1.27, the usual axial pion-nucleon coupling. 
Note that in the MSSM, c v = c a = 0. 

N 2 X 2 vertex. The one-nucleon contribution to the 
scattering amplitude is traditionally the only term con- 
sidered and the full N 2 X 2 Lagrangian is given in Ref. 
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In contrast to the two-nucleon case, the sum in Eq. (JTJ 
is now over all quark flavors. At LO and ignoring the 
pion-pole term for simplicity, we have 



£ iV £ SN(d [ > + d L s T d )N + V^N^idl + dy)N 

+A^l\dl + dlr^+T^a^ + dlr 3 ^, (10) 

where [with (k,l,m) — (1,0,0)] 



where M is the nucleon mass, S^nnH^) — 13.7 is the 
strong pion-nucleon coupling constant and the nuclear 
operator is given by 

Oq v {xx, ...,Xi) = -5(x! - x 3 )S(x 2 - xa) 
x-[F 1 (x)a a p ■ a^s + ¥ 2 {x)T a<j) ^ s \f ae -f 7C , (18) 



4 2 



-2$> 

s,b,c,t 



d i=— («"- 



^4,(11) 



and j — {1, 3, 6, 7} for i — {s,v,a,t} respectively. The e's 
are 0(1) except for el = (N\uu—dd\N) — where N can be 
either a neutron or a proton — which is small numerically 
and will henceforth be neglected. For the scalar terms, 
we have in particular [ill Il2| 



e° = (N\uu + dd\N) = 16 ±8, 
e s s = (N\ss\N) « (0.04-0.2)e s °, 



2M 



m d 



M 



(12) 
(13) 

,(14) 



where Q = c,b, t. Although Ref. 12] suggests the larger 
value for e| ~ 0.2e|?, the error is very substantial and we 
will use the more central value e| ~ O.le". 

MSSM. In the MSSM, the WIMP is the Majorana 
spin-1/2 neutralino which couples to quarks according to 
the low-energy Lagrangian of Eq. Q with a| = a\ — 
a 5 = a 7 ~ |13| . Since we are interested in compar- 
ing the one- and two-nucleon scalar-scalar interactions of 
Eq. we will not consider the axial-axial a\ term and 
will only focus on a\ which in the MSSM is given by 



M w 



(15) 



where m q is the quark mass and the Sq's are dimen- 
sionless and are usually larger than 10 in the MSSM [l3| . 
The LO SI x 2 -hadron Lagrangian generated by the scalar 
currents is simply 



Xhad 

— 



(d° s NN + b s lf ■ 7?) XX- 



(16) 



Therefore, in the MSSM, an important modification of 
the SI x 2 -hadron Lagrangian is the new 7T 2 x 2 term at LO. 
This term, through the two-nucleon diagram of Fig. [Ja), 
can give a contribution comparable in size to the one- 
nucleon scattering amplitude. 

The calculation of the diagram in Fig.^a) is straight- 
forward 0: we use the Feynman rules to evaluate the 
amplitude for two particular nucleons inside the nucleus 
and we take the non-relativistic limit; we then Fourier 
transform the result and obtain the operator 



. 5*7,- ran g^jyjv 
l 3V2M 2 4tt 



where the £'s and r/'s are spinors and isospinors respec- 
tively, T^ps = ia a p ■ pa^s ■ P - v a p ■ o^a and 

Pi(x) = (x-2)e- x , F 2 {x) = (x+l)e- x , (19) 
S, = ( S u ^ + S d ^-) , (20) 



(m u + m d ) 



w 



Mi 



w 



where x — m^p is proportional to the distance between 
the two nucleons with p — x\ — x 2 , p — \p\ and p = pj p. 
In Eq. (|18J) . a and f are spin and isospin Pauli matrices 
respectively. 

The contribution to the WIMP-nucleus scattering will 
therefore be 



A" 



M 
1 x 



3V2M 2 4tt 



(Z,A\-^0^\Z,A), 



(21) 
(22) 



where the sum is over all nucleon pairs inside the nucleus, 
A is the number of nucleons, and Z is the charge of the 
nucleus. 

On the other hand, the one-nucleon contribution will 
simply be 



A NN = 2id°A. 



(23) 



Taking the ratio of the two-nucleon to one-nucleon con- 
tribution we obtain 



~(0.21±Q.08)r^,r-^ + ^ 



A NN 



A 



GjNN 



(24) 



S NN = \{S u m u + S d m d )e° s + £ S, 



q m q el (25) 



■jj.t 



Note that we did not include contributions from heavy 
quarks with extra, inverse powers of the squark masses 
for simplicity , and that the error in the overall factor 
of 0.21 comes from the uncertainty of the current-quark 
masses. Looking at the special case where all the S q are 
equal and cancel in the ratio r with m s = 115 ± 55 MeV, 
we find 0.002 < (0.21 ± 0.08)r < 0.006. Since there 
exists regions of the MSSM parameter space where the 
S q s have opposite signs, the ratio r can be greater than 
one if cancellations occur. For example, with S u = S d = 
S s = -S c = -S b = S t , m s ef = 115 MeV, and (m u + 
m d )/2e° = 60 MeV, r = 1.5. Alternatively, if all the S q 's 
have the same sign, a cancellation could still occur with 
the extra heavy quark terms not included in Eq. (|25|l if 
they are of similar magnitudes but opposite in sign. 
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We estimated AT™ in the independent particle ap- 
proximation assuming closed shells and N = Z where 
N is the number of neutrons. We performed the calcu- 
lation in two different mean- field potentials: a square 
well and a harmonic-oscillator potential. To account 
for the hard core, we used Jastrow functions {defined 
[1 — exp(— ap 2 )(l — bp 2 )] with a = 1.1 fm~ 2 and b = 
0.68 fm~ 2 } as well as a simple integral lower cut-off at 
0.5 fm. We found that the matrix element is approxi- 
mately linear in A for nucleons in a harmonic potential 
and exactly linear in A for nucleons trapped in a square 
well. We also found a minimum value of Af™ — A in 
the harmonic oscillator potential with Jastrow functions 
and a maximum value of A/q™ = 2A for the square-well 
potential with the lower cut-off at 0.5 fm. We estimate 
that ATq^ I A w 1.5 ±0.5 although the error can be signifi- 
cantly larger due to the fact that the nuclei used in direct 
dark-matter searches are not closed-shell nuclei and that 
they have a neutron excess. We also do not take into 
account spin-orbit coupling. Furthermore, Affi^/A could 
vary substantially from one nucleus to the next. Keeping 
these issues in mind, we conclude that Af™ and A are of 
the same order and that the ratio in Eq. l)24f) can be of 
the order of one or larger depending on the values of the 
MSSM parameters. 

In supersymmetric models, the one-nucleon current 
generically produces roughly equal SI couplings to the 
proton and neutron [5j, which results in a SI amplitude 
that is proportional to the atomic number of the nu- 
cleus. Inclusion of the two-nucleon contributions could 
change this picture since such contributions might can- 
cel against the one-nucleon contributions. If the ratio of 



the two-nucleon matrix element to the atomic number 
varies from one nucleus to the next so will the degree of 
the cancellation. Thus, when the two-nucleon contribu- 
tion is taken into account, a dark-matter candidate that 
appears in DAMA but not in other searches 0] is con- 
ceivable for a WIMP with SI interactions even within the 
framework of the MSSM (Ref. [r| provides arguments 
against a WIMP with SD interactions as an explanation 
for the DAMA signal). 

With 1 < Af™/A < 2 and r m 1, the exchange- 
current contribution is in the range of 10-60%. This 
is atypical since such currents usually contribute about 
5-10% to the amplitude asinp + n— » ci + 7 and 
vi+d~^v+p + n 17] . Part of the reason the exchange 
current contribution is normally small at low-energies 
stems from the fact that the "probe" -tt 2 vertex (where 
the probe can be a photon or a neutrino for example) is 
typically suppressed by at least a factor of q/M ~ 0.1. In 
WIMP-nucleus scattering, the situation is different since 
the tt 2 X 2 vertex need not be suppressed by this factor as 
we showed in Eq. (|16|l . Thus, instead of a 5-10% effect on 
the WIMP-nucleus scattering amplitude, the one-nucleon 
and two-nucleon contributions could be comparable. 

In conclusion, we have shown that the two-nucleon and 
one-nucleon currents can give comparable contributions 
to the WIMP-nucleus scattering and point out that this 
result can help resolve the conflict between the various 
direct dark-matter search experiments. 
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9821 and DoE DE-FG03-92-ER40701 and DE-FG03- 
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